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Abstract 

We study the generalized Chaplygin gas (GCG) scenario in a modified 
gravity approach. That is, we impose that our universe has a pure dust 
configuration, and allow for a modification of gravity that yields a GCG 
specific scale factor evolution. Moreover, assuming that this new hypothetical 
gravity theory obeys a generalization of Birkhoff's law, we determine the 
Schwarzschild-like metric in this new modified gravity. We also study the 
large scale structure formation in this model. Both the linear and non-linear 
growth are studied together with the growth of the velocity fiuctuation in the 
linear perturbation theory. We compare our results with those corresponding 
to the ACDM model and discuss possible distinguishable features. 
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1 Introduction 



The nature of the dark energy (DE) that is responsible for today's cosmic 
acceleration is an open and tantalizing mystery ^J. It leaves room for novel 
theoretical explanations and new cosmological scenarios. The most obvious 
is the A Cold Dark Matter model (ACDM), comprising dark matter and 
a non zero cosmological constant-A. But the difficulty of explaining the 
value of this A term from fundamental physics leaves motivations for other 
phenomenological proposals. Two such alternative avenues for DE model 
building in which the DE is dynamical are the quintessence [2] and k-essence 
[S] models. Similar to inflation, the first one is constructed with an ordinary 
minimally coupled scalar field whose equation of state is a function of time. 
Amongst them, the tracker quintessence models j4j have the advantage of 
allowing the current accelerating epoch to be reached from a large set of 
initial conditions. 

Recently, there are suggestions that the present acceleration of the uni- 
verse is not due to any new unknown component in the cosmic soup, but 
due to a modification of the gravitational physics at scales typically much 
smaller than today's horizon. A few suggestions are particularly interesting 
in this regard: the Dvali-Gabadadze-Poratti (DGP) brane induced gravity 
model the Cardassian model proposed by Freese and Lewis jHI and also 
the recent model by Dvali and Turner [7j. All these models lead to late-time 
acceleration with a modified Friedman equation and with no explicit dark 
energy component. The basic feature of this approach is that the geodesies 
in a static spherically symmetric spacetime may completely determine the 
cosmological evolution, i.e., the present acceleration of the universe is not 
driven by any extra dark energy component but by the matter content itself. 
Given a complete cosmological evolution of our universe as suggested by ob- 
serv! ations (i.e. given a specified scale factor a(t)), if we impose that it is due 
to a pure dust configuration, we can study what modification of the gravity 
is necessary to yield such a scale factor. Moreover, imposing that this new 
hypothetical gravity theory obeys a generalization of Birkhoff's law, one can 
get a unique modification. This procedure of determining the Schwarzschild- 
like metric of the new modified gravity that gives the prescribed cosmological 
evolution with a dust-filled universe has been studied extensively by Lue et.al 
jH]. Also, in a recent paper Multamaki et.al [§] have studied the growth of 
large scale structure in the DGP and Cardassian model. They put forward 
a general formalism for calculating the growth of both linear and non-linear 
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fluctuations in models with non-standard Friedman's equations. 

Another alternative to the quintessence model which has attracted great 
interest in recent times is the so-called generalized Chaplygin gas (GCG) 
model [ini- The model explains the acceleration of the universe via an exotic 
equation of state resulting in a behaviour like dark matter at early times and 
dark energy at late times. The GCG is characterized by the equation of state 



where A and a are positive constant. For a = 1, the equation of state is 
reduced to so-called Chaplygin gas scenario [TT] . 

Inserting the above equation of state in the energy conservation equation, 
one can integrate it to obtain 



where B is an arbitrary constant of integration which should also be positive. 
One can see at once that this energy density interpolates between a dust like 
configuration in the past and a de-Sitter like one in the late times. This 
property makes the GCG model an interesting candidate for the unification 
of dark matter and dark energy. 

Using the above expression for p in the Einstein equation for H{t), one 
gets 



This expression for H{t) is a very good fit for the observational data as far 
the background cosmology is concerned. This has been successfully con- 
fronted with various observational test: high precision Cosmic Microwave 
Background data jT2], Supernova data 13^ and data from gravitational lens- 
ing [Hj. 

But despite all these pleasing features, the main difficulty with such an 
unified model is that it produces unphysical oscillations or exponential blow- 
up in the matter power spectrum at present ^3] . Moreover, it was shown that 
the linear approximation breaks down at an early stage, implying that a more 
careful approach, including nonlinear effects, should be taken into account 
Some efforts have been made to circumvent the previous problem by 




(2) 




(3) 
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adding the baryons into the model which is not accounted for by GCG pTj . 
GCG has also been treated as a dark energy model in combination with 
dark matter, where only the dark matter part is perturbed ^H]- In a very 
recent work, it has been shown that this equation of state represents uniquely 
an interacting mixture of decaying dark matter and a cosmological constant 
once one excludes the possibility of having a phantom like dark energy and 
in such a scenario one can avoid the problem of having unphysical features 
in the matter power spectrum [T^. 

In this letter, we are considering the GCG in a modified gravity approach 
as described by Lue et.al [S]. That is, we assume that the background evo- 
lution is due to some kind of modified gravity, rather than an energy density 
with a GCG equation of state. We impose the same background evolution, 
but with an energy density consisting purely of matter (we ignore the residual 
radiation term at present). We deduce the corresponding Schwarzschild-like 
metric for this new modified gravity. We also study both the linear and sec- 
ond order density perturbations in this modified GCG universe and compare 
our result with that of a standard ACDM model. 



We assume that our background universe is well described by the H{t) given 
by Eq. Elbut that it only contains a pure matter configuration. That is, our 
assumption is that the only energy density in the universe has an equation 
of state w = with a conservation equation 



Subsequently we look for the possible modification of Einstein gravity which 
results in the above H{t). Noting that for the dust in Eq. HI p(t) oc a~^{t), 
one can see that such modification is given by 



One can relate the constant A with the matter density parameter fi^ as A = 
(3i/oV87rG)("+^)(l Also, we have identified the previous constant B 

with Pq~^^ (po being the present matter energy density), which can be done 
with no loss of generality. Notice that since we asume our energy density 
restricted to dust, we have Qt = < 1- Alternatively, we can interpret 



2 Modified GCG model 



p = -3Hp 



(4) 
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the constant A as a cosmological constant, adding an extra component to the 
energy density, hke what is done in Either way, the evolution equation for 
the matter density Eq. Hlis unchanged and does not couple to any other fluid. 
This is in contrast to what was done in ^H]- (Furthermore, when a = this 
expression reduces to a ACDM model. Since we are only considering matter 
as our energy density, this means that the two models will have exactly the 
same evolution. That is, for a = we have a "modified gravity" ACDM 
model). ^ Using Eq. one can write the modified Einstein equation as 

H^ = Hlg{x), (6) 

where g{x) = [{1 — Q'^^) + x°'^^f'^"'^^ and x is a dimensionless quantity 
defined as x = p/Hq. For x >> 1, g{x) x and one recovers the 
standard Einstein equation in the early universe, whereas in the late time 
the gravity is modified. Since we are assuming the scale factor evolution to 
be due to a modification of Einstein's gravity, it is interesting to consider 
what effects this modification could have at astrophysical scales. This is 
particularly easy to do if we impose an extra condition on the theory, namely 
that it obeys a generalization of Birkhoff 's law. Following the procedure in 
jHl, it is possible to deduce the Schwarzschild-hke metric of this modified 
gravity around a spherically symmetric matter source. In general, this is 
given by 

3 

goo = 9;r=l-r'H^9C-^), (7) 

where g is the function defined after Eq. r is the usual radial coordinate, 
and is a measure of the distance scale over which gravity is modified. It 
is given by = {^^^Y^^, M being the mass of the spherically symmetric 
gravitating object of energy density p and radius r and is given by M = 
^p{t)r^. Since g[x) x when x ^ 1, and we have x = ^, it easy to check 
that the metric will go to the usual Schwarzchild metric at small r -C r^. 
However, for values of r ^ Tc the metric will be modified, and one then 
expects to see a significant deviation from the usual Einstein gravity. As an 
example, consider the universe within the horizon, with a mass M -(T^- 

■'For the usual GCG scenario, the a = case is also equivalent to a ACDM model, 
however one has to be careful with the background used to calculate the perturbations 
when comparing between the two j20|. 
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This will give us a scale ^ Hq , as expected, since it is approximately the 
scale at which acceleration sets in. 

Even if we consider radii smaller than Tc, one can still try to use this 
metric, with its small modifications, to impose some astrophysical bounds 
on this model. As a particular example, for the solar system the orbits of 
the planets are well within the value of Vc- However, we can expand the 
metric in powers of — and check what the first order correction to Einstein's 
gravity gives. This was done in , where the authors obtained a formula for 
the precession of the perihelion. For our modified gravity, the correction will 
be proportional to (r/rc) ^^^^^^^^ If the exponent is negative, the correction 
would be huge for very small values of r. This imposes a bound on our 
parameter space of a > — |. In our numerical simulations we will restrict 
ourselves to a > 0. 



3 The perturbed equation 

To study the density perturbations, we consider the evolution of a perturbed 
ideal fiuid, with a shear free four- velocity u'^. We choose a coordinate system 
such that u'^ = (1, ax + v), with v being the peculiar velocity of the fiuid. 
We define the perturbed part of the energy density 6 as 

p = p(l + 5) (8) 

where p is the background energy density (and for the rest of this paper 
barred quantities refer to their background values). The perturbed part of 
the continuity equation is then given by j|21j 

^ + (1 + 6)6 = (9) 
ar 

where r is the conformal time, dt = adr, and 6' = V • v. 

We can get a second equation for the fiuid evolution from Raychaudhuri's 
equation. This can be written as [0] 

6 9 

- + -H = 3{H + H"^) -3{H + H^) . (10) 
Combining the two, we get an equation for the perturbed energy density 
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Figure 1: Linear growth Di for a = 0, 0.2, 0.4, . . . , 1.0 from bottom to top. 



evolution, for a general matter density background, 

(11) 

Following [Oj, we perform an expansion in the perturbation 5 

°° °° D (n) 
5 = Y^5, = Y^^5l (12) 

i=i i=i 

where 6q is a series expansion parameter. We will be interested in the first 
two terms in the expansion, for which the linearized equations are, in general, 

D'i + (2 + ^)D[ + SciZ^i = (13) 

^2 + (2 + §^)D'^ - \{D[Y + 3ciD2 + 6(ci + c^)Dl = (14) 

It is straightforward to compute the two terms Ci and C2 for the "gener- 
alized Chaplygin" scenario. They come out as 
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~ 2(A + p"+i)2 

^2 = 4(;+',l\)3 [-a + 3")-4+(2 + 3c.)p""^] (16) 

and it is easy to check that they reduce to the standard cosmological constant 
result for a = 0, as expected. The second order perturbation is related with 
the skewness of the density field at large scale. The q-order moment of the 
fluctuating field is related with the perturbation as 

ruq =< 5^ > . (17) 

The normalized skewness is given by 

S3 = ^ (18) 

1712 

which can be written in terms of the first and second order perturbations. 
For Gaussian perturbations < 5^ >= 0, so that one gets 

^3 = 3^ (19) 

In the standard CDM model, this coefficient can be calculated exactly to 
give ^3 = 34/7 ^ 4.86. 



Numerical results 

At early times, there is no difference in the evolution between the standard 
CDM and any of the proposed modified gravity scenarios, including ours. 
Therefore, for our numerical results, we start the simulation at an early time, 
a = 10~^, taking as initial conditions the standard CDM solution, S oc a. In 
calculating the second order perturbation, initial conditions are chosen such 
that the standard solution 5*3 = 34/7 is valid from the beginning. In figure 1, 
we have plotted the first order perturbation Di for different values of a. In 
Figure 2, we have shown that skewness S3 for the second order perturbation. 
In plotting these figures, we have assumed a fixed = 0.3, which sets the 
value of A for each a (cf. the comment after Eq. [SJ. 
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Figure 2: Non-linear growth 5*3 for a = 0, 0.2, 0.4, . . . , 1.0 from top to bottom 
(at present). 
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Figure 3: Variation of the ratio /mg//a for a = 0, 0.2, 0.4, . . . , 1.0 from 
bottom to top (at present). 
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Figure 4: Contours for the present day value of / in a-^lm plane. 



The general form of linear growth Di is similar to the ACDM case but one 
can have more linear structure formed at present with increasing values of 
a. This will also result in a higher ag value (the rms fluctuation on a sphere 
of 8 Mpc/h) at present than the corresponding ACDM value. Also one does 
not have any oscillations or exponential blow in the matter power spectrum 
at present as in the usual GCG scenarios. This is expected as a consequence 
of having only a dust configuration in the universe with a vanishing pressure. 

In contrast, the second order perturbation or 5*3 evolves in a completely 
different way from the ACDM model, as shown in figure 2. However the vari- 
ations from the ACDM case at present are approximately 2% for a as high as 
1. Current estimates for agree with the standard prediction but with large 
uncertainties, of the order of 20% — 30% [22j. Hence the current observational 
result can not be used to differentiate between this modified GCG gravity 
with the standard ACDM model using S^. Although the Sloan Digital Sky 
Survey (SDSS) is expected to reduce these uncertainties to around 5%, still 
this will not be enough [231 • have also checked that varying Qm within 
the range 0.2 < fi^ < 0.4 does not change the result much. 

In figure 3, we have shown the variation of the ratio /mg/Za, where / = 
d InDi/d In a. One can see from figure 3 that for a = 1, f can increase 
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up to 40% from the corresponding ACDM value. Notice that / governs 
the growth of the velocity fluctuations in the linear perturbation theory; 
therefore the large deviations of / with changing a is detectable via precision 
measurements of large scale structure through joint measurements of the 
redshift-space power spectrum anisotropy and bi-spectrum from z = to 
z ^1. The SDSS should be able to probe this quantity with statistical error 
of order of a few percent j2Sl- In figure 4 we have also shown the contour plot 
for the present value of / in the a-f^m plane. It shows that with increasing 
fim, the / depends more strongly on a. 

4 Conclusion 

Using as motivation the good observational fit of a GCG model to the back- 
ground evolution of the universe, we have studied it in a modified gravity 
approach as described by Lue et.al |[^. We have also assumed that the new 
gravitational physics obeys a generalization of the Birkhoff 's law. This means 
that an observer in the gravitational field of a spherically symmetric source 
of mass M experiences a significant deviation from the usual Schwarzschild 
metric at a distance scale greater than approximately (^^7^)^''^, where i^o is 
the present Hubble radius. /^From this, a simple calculation with the metric 
Eq.[7| shows that astrophysical bounds impose a value oia> —\. An inter- 
esting check to the validity of this model would entail a more detailed study 
of its behaviour at astrophysical scales. 

We have studied the first and second order density perturbation in this 
model. Our results show that the linear perturbation Dx evolves in a similar 
way to a ACDM model, but with a larger value at present. This results 
in an enhancement in the corresponding cxs value, namely up to 22% at 
present for a = 1. In theory, the nonlinear perturbation 6*3 gives a more 
radical signature, evolving quite differently from the ACDM case. However, 
the effect is very small, and with the present observational uncertainty, one 
can still not use this to distinguish between the two models. We have also 
studied the parameter / related with the velocity fluctuations in the linear 
perturbation, and have shown that the variation obtained in comparison to 
a ACDM model is quite signiflcant. This leads to a possible detection by the 
present day observations. 

As expected, by considering a modifled gravity approach to a GCG, we 
can avoid the problem of having unwanted oscillations or exponential blow- 
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up in the matter power spectrum at present, as one expects in its unified 
dark matter-dark energy approach. 

Finally, it remains to be seen how one can ultimately obtain such modi- 
fications in gravity from a fundamental theory. But given the fact that the 
GCG type equation of state can arise from Born-Infeld type lagrangian (TU], 
one may expect that D-brane physics can shed possible light to tackle this 
problem. 
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